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Abstract — An unsteaih two dimensional magnetohydrodynamic flow of a visco-elastic incompressible fluid 
(Walters B' fluid model) along an infinite hot vertical porous surface hounded by porous medium with free 
stream as well as suction velocity in the presence of heat surface and chemical reaction lias been investigated. 
The governing equations of motion, energy and concentration are solved by the successive perturbation 
technique. The flow phenomenon has been characterized by magnetic parameterM , elastic parameter K Q , 
porosity parameter K p , Grashoff number for heat transfer G r , Grashoff number for mass transfer G c , Prandtl 
number P r and Schmidt number S c . The expression for skin friction coefficient, Nusseli number and Sherwood 
number at the surface are also derived. The variations in the fluid velocity, temperature and concentration are 
shown graphically where as numerical values of skin-friction; Nusselt number and Sherwood number are 
presented in tabular form. 

Key words — Chemical reaction. Heat source, MHD flow. Porous medium, Walters 's B' model. 

I. Introduction 

The influence of magnetohydrodynamics heat and mass transfer processes over a moving surface is of 
interest in engineering and geophysical applications such as geothermal reservoirs, thermal insulation, enhanced 
oil recovery, packed-bed catalytic reaction and cooling of nuclear reactors. Many chemical engineering 
processes like metallurgical and polymer extrusion processes involve cooling of a molten liquid being stretched 
into a cooling system; the fluid mechanical properties of the penultimate product depend mainly on the cooling 
liquid used and the rate of stretching. Some polymer fluids like polyethylene oxide and polyisobutylene solution 
in cetane, having better electromagnetic properties are normally used as cooling liquid as their flow can be 
regulated by external magnetic fields in order to improve the quality of the final product. 

Flows through porous media are frequently used in filtering of gasses, liquid and drying of bulk 
materials. In electrochemical engineering, porous electrodes and permeable, semi permeable diaphragms are 
used to obtain improved current efficiencies. In the field of agricultural engineering, porous media heat transfer 
plays an important role particularly in germinations of seeds. Above all man does a part of his breathing through 
his porous skin. More interestingly, application of magnetic field in the field of chemical engineering and 
agriculture attract many researchers to work on this new field of interest. Kronenberg rightly stated that "if the 
agricultural industry knew the benefit of magnetically treated water, the total agricultural output could be 
approximately increased by 25%". When fluid passes through a properly focus magnetic field, the magnetic 
field breaks up some of the complexes that are carried in the water, freeing the captive mineral particles. This 
process keeps minerals in the water rather than precipitated out. This is known as magnetic field conditioning. 
Further, magnetically treated fuel has a tendency to attract oxygen molecules when mixed with air in a 
combustion cylinder. This results a more efficient and complete combustion of the fuel generating more power 
from same amount of fuel. Research has shown that dairy cows will drink more magnetically treated water that 
untreated. This has increased milk production appreciably. With beef cattle, feed conversion ratio has increased 
by 5% approximately enhancing meat production about the same. A chemical reaction involves in breaking of 
bonds in the reactive substances and formulation of bonds to form product species. 

Several researchers have studied the two dimensional free convection, heat and mass transfer flow of 
an elastico- viscous fluid through porous medium. Sharma and Mathur [1] have studied the steady laminar free 
convection flow of an electrically conducting fluid along a porous hot vertical plate in the presence of heat 
source/ sink. Dash and Rath [2] have studied the problem of laminar flow and heat transfer of an electrically 
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conducting fluid between parallel porous plates by applying explicit finite difference scheme. Walters [3] has 
studied the Non-Newtonian effects in some elastico-viscous liquids whose behavior at small rates of shear is 
characterized by a general linear equation of state. Unsteady two dimensional laminar flow of a viscous 
incompressible electrically conducting fluid in the vicinity of a semi-infinite vertical porous moving plate was 
discussed by Kim [4]. Flow and heat transfer of an electrically conducting visco-elastic fluid between two 
horizontal squeezing/ stretching plates has been studied by Rath et al. [5]. Unsteady MHD free convection and 
mass transfer flow past an infinite heated porous vertical plate with time dependent suction was discussed by 
Cookey and Sigalo [6]. 

Ece [7] has studied the free convection flow about a cone under mixed thermal boundary conditions 
and a magnetic field. Pravazinia and Nassehi [8] have studied the sheer thinning fluid flow through highly 
permeable porous media. Free convection heat and mass transfer from a horizontal cylinder of elliptic cross 
section in micro polar fluids has been studied by Yang [9]. Unsteady flow and heat transfer through an elastico- 
viscous liquid along an infinite hot vertical porous moving plate with variable free stream suction have been 
studied by Sharma and Pareek [10]. MHD flow through a porous medium past a stretched vertical permeable 
surface in the presence of heat source/ sink and a chemical reaction have been studied by Dash et al. [11]. The 
effect of Hall current and chemical reaction on MHD flow along an exponentially porous flat plate with internal 
heat absorption/generation was studied by Rath et al. [12]. Free convective MHD flow through porous media of 
a rotating oldroyd fluid past an infinite vertical porous plate with heat and mass transfer have studied by Dash 
et al. [13]. Sharma and Sharma [14] have studied the unsteady two-dimensional flow and heat transfer through 
an elastico-viscous liquid along an infinite hot vertical porous surface bounded by porous medium. 

The main objective of the present problem is to investigate the combined effect of transverse magnetic 
field, heat source and chemical reaction in the presence of oscillatory suction and free stream velocity on the 
flow, heat and mass transfer phenomena of an electrically conducting visco-elastic fluid. 

II. Formulation of the problem: 

An unsteady two dimensional MHD flow past a vertical infinite surface through a porous medium of a 
visco-elastic fluid (Walters fluid (model B' )) has been considered. The time dependent fluctuating suction 
velocity has been introduced with an oscillatory free stream velocity. A uniform transverse magnetic field B is 
applied normal to the direction of the fluid flow. The X - axis is taken along the surface in upward direction i.e. 
opposite to the direction of gravity and y -axis is taken normal to the surface. 

The equations of continuity, motion, energy and diffusion for flow of a visco-elastic fluid through 
porous medium bounded by an infinite, hot vertical porous surface with oscillatory suction velocity are given by 
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Where U , V are the velocity component along X and y direction respectively , p the density of 
the fluid , g the acceleration due to gravity, ju the dynamic viscosity, K p the permeability parameter, 
K the non-Newtonian parameter, K the thermal conductivity of the fluid, C P the specific heat of fluid at 
t pressure , K c the chemical reaction parameter, S * the heat source parameter , T the temperature 
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of the fluid, J3 and J3 are the volumetric coefficient of thermal and concentration expansion respectively, a 

the electrical conductivity , C the concentration of the fluid , D the mass diffusion co-efficient and t the 

time. 

The boundary conditions are: 

y * =Q:u* =Q,T* =T w ,C* =C W (5) 

y* ->oo: M * -^U*(t*]T* -^T X ,C* -*C X 

Where T w is the surface temperature, U*[t*j the free stream velocity, T x is the free stream temperature, C w 
is the concentration at the wall and C x is the concentration of the fluid far away from the wall. 

From equation of continuity (1) it is clear that the suction velocity normal to the plate V is either a 
constant or a function of time. Hence, it is assumed that V* is in the form of 

v* = -v {) (\ + se ic ° r ) (6) 

Where CO is the frequency of vibrations, e is the small parameter i.e., (0<£ <1) and V is a non zero positive 
suction velocity. The negative sign indicates that the suction is towards the plate. 

To reduce the governing equations into non dimensional form, let us consider the following 
transformations. 

= /v t j\ 2 u = u a= 4va U(t) = U'(f) T = T'~T X p = juC P R = K'v 2 
pv v D v Q y vp C w - C x v 

Gm «rfr.-T.) ^-(c.-c.) 

Uv Q ' Uvl 

Where G r , G c , P r , S c , K p , K c , M and K Q are the Grashoff number for heat transfer, the 
Grashoff number for mass transfer, the Prandlt number, Schmidt number, the porosity parameter, the chemical 
reaction parameter, the magnetic parameter and the non-dimensional non-Newtonian parameter. 

Using the transformation (7) equation (6) reduces to: 

v'=-v (l + *<?**) (8) 

where CO is the frequency of the suction velocity. 

In view of equation (7), equation (2) - (4) become 

1 dw ,. im.du d 2 u .. icot . „ \\ d 3 u ,. i01t .d 3 u} 
(l + ee ) — = — v + m(l + £e -u)-KA T -(l + ee ) — ^ 
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where m = M H (say) . 

K p 

The corresponding boundary conditions in non- dimensional form are 

y = 0:u=0,T = l,C = l 

y -> oo :u ->[/(/), :r->0,C->0 (12) 

In the neighborhood of the surface it is assumed that 
u(y,t)= u Q (y) + £u 1 (y]e" n ' , 

T{y,t) = T (y) + sT l {y)e iat , (13) 

C(y,t)=C (y) + eC l (y)e ia ' 
and for the free stream U(t) = 1 + «>"* . 

Now, substituting (13) in the equations (9)-(ll) and equating the coefficients of like power of £ we 
get the following zeroth- order equations: 

K q u"' + u'q +u' q —mu = —G r T Q —G C C —m (14) 

T;'+PX + SP r T Q =0 (15) 

C" + S C C' Q - K C S C C = (16) 

First order equations are: 

K u 1+ (l - ^ V + u[ -L + ^ V = -G r T x - G C C X -[m + ^y K u m - u' (17) 

77+ p r r; + [ p r s - ^- V - -p r r ' ( 1 8) 

C"+S c C[-\ -^— + K C S C C, --5 C C; (19) 

where prime denotes the differentiation with respect to y . 
The corresponding boundary conditions are changed to 
y = 0:u =0,m, =0,T Q =1,7:, =0,C =1,C, =0 

j -> oo : w -» 1, M] -> l,r -> 0,7; -> 0,C -> 0,Q -> (20) 

Solving equations (15), (16), (18) and (19) using the boundary conditions (20), we get 

T = e m ' y , (21) 

C = e m " y , (22) 

T^I^-e^) , (23) 
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C 1 = I 2 (e m " y -e m " y ) (24) 

Pm 2 4im 6 
where l y = and I 2 = 

r^+P T m 2 +SP r - — P T ' a 

2,2 r 4 , 

But equation (14) and (17) are of third order equations when K Q ^ and these equations reduce to 
second order equations when K = (for Newtonian fluid). Hence the presence of non-Newtonian parameter 
increases the order of the differential equation. As the non-Newtonian parameter (K Q ) is very small for 
incompressible fluid (Walters), therefore U Q and W, can be expanded in powers of K Q which are given by 

u {y) = u 00 (y) + K u 01 (y)+o(K^) 
u 1 {y)=u w {y) + K u ll {y)+o{K 2 ) 

Introducing (25) into the equations ( 14) and (17) and equating like powers of K Q in both sides we get 

M oo + M oo _ mM oo = ~GT a - G C C Q - m (26) 

«,"n + u'-\m + — L in = -G T, -G r C, -\ m + —\- u' (27) 

K\ + M oi _ "^oi = _ "oo ( 28 ) 

M - +u[i J m + i^\ =_ <o +<o !^_ M » - Ki (29) 

ii ii ^ 4 j u io io 4 00 01 

The boundary conditions are 
y = : u 00 = 0,w 01 = 0,u w = 0,u n = 

y -»oo:w 00 -»l,w 0] -»0,w 10 -»l,w n -»0 (30) 

Solving equations (26)- (29) using the boundary condition (30), we get 

w 00 = I 3 e m '" y - I 4 e miy - I 5 e m °' + 1 (31) 

u 01 = I 9 e m '° y + I 6 ye mioy + I 7 e m2y + I 8 e m6y (32) 

"io = h^ mny + hi^ ~ V" 4 " + h^ ~ h^ y ~ hie miay + 1 (33) 

M n = h5 e " h2y + hs em2y + ^29 e, " 4 ' V + ho^^ + h^""^ + hi 6 " 1 ^ + h^'^y ~ h^™'^ y (34) 

Hence the velocity w(y,?) is given by 

w(;y,0 = I 3 e m '" y -I A e miy -I s e m " y +l + ^ {/ 9 e m|0> ' +I 6 ye m, " y + I 1 e miy + I 8 e m<,y \ 
+ s\l,,e" h2y +L 1 e miy -Le'" t,, + I w e mt,y -I on e mi>y -I 2] e mi ° y +1 

L16 (35) 

+ K Q [I 35 e m ' 2y +I 2 ^e" hy +I 29 e'" i} + I 3Q e mf,y + I 31 e m * y + I 32 e m '° y 

+ I 33 e mny y-I 34 e mioy y}y wt 

Where I 3 —I 35 ,tn l , ,m 12 are given in the appendix. 

Skin friction [C f J 
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c f =-^ = r*i-jr (l^-(i + «'-)^Y| 06) 

/ pC/u [ej \4 9% V ' By J_| =0 

Nusselt Number [N u ) 

N u = «» , = -^1 =-[m 2+ ^/ l( m 4 -m 2 )] ( 37 ) 

" v K(T m -T°o) By \ y=0 L 2 lV 4 2 J 

Sherwood Number (S A ) 

5/1 V(|jD (c„ - c, ) ay J y=0 Lmf ' + * /2 (m ° m * )J 

Where T ffl is the dimension shear Stress component of the elastico-viscous fluid and the symbols q x and 
q 2 represent heat and mass flux and these are given by 



l^dy* 



III. Results and Discussion: 

The following discussions bring out the effects of porosity parameter(A" P ), Source parameter(S), 
Grashoff number for heat transfer ((>',), Grashoff number for mass transfer (G c ), prandatl number(fV), Schmidt 
number(,SV), chemical reaction parameter(/Q, magnetic parameter(M). phase an«le(cjf) and elastic parameter^). 

Fig. 1 shows the velocity distribution for various values of G„ G a K , K p , P„ S c , co, cot, and M. The 
common characteristics of all the curves are to increase the velocity with the layers near the vertical surfaces. 
From the curves (III, IV, V, VI) it is observed that an increase in G, increases the velocity near the plate. It is 
due to free convection current which accelerates the velocity near the plate. Further higher Prandtl number fluid 
contributes very insignificantly in decreasing the velocity (curves III, VII). It is further to note that the presence 
of transverse magnetic field reduces the velocity due to generated Lorentz's force opposing the motion (curves 
III, XIII). It is also observed that fluid velocity increases due to increase in the elastic parameter(A" ), Grashoff 
number for mass transfer (G c ), chemical reaction parameter(A^ c ), and Source parameter(S) but reverse effect is 
observed in case of porosity parameter(A"p), phase angle(»0» an d Schmidt number(5 c ). 
Curve GPKCJwt KGMKS S 
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Fig- 1 Effect of different parameters on Velocity profile when s = 0.5 

Fig. 2 depicts the temperature distribution in the flow domain. An increase in heat source parameter (S) 
contribute to increase the thickness of the thermal boundary layer (curve I, V, VI) but with an increase in the 
value of prandtl number (P r ), frequency parameter (co) and phase angle (cot ) the temperature decreases (curves 
I, II, III and IV). This may be attributed to the fact that smaller values of P r lead to increase the thermal 
conductivity of the fluid and there by more amount of heat is diffused which causes the decrease in temperature. 
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Fig- 2 Temperature profile when U = 0.5 

Fig. 3 exhibits the concentration distribution in the flow field. An increase in Schmidt number (S c ) and 
chemical reaction parameter (K c ) leads to decrease the concentration. Also an increase in frequency of 
oscillation^ leads to increase the concentration (curves I, VI). It is also interesting to note that phase angle has 
no impact on the concentration profile as all the curves coincide with each other (curves I, VII). 
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Fig- 3 Concentration profile when U =0.5 

Table -1 Skin friction (C f ) 
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Table 1 shows the numerical values of skin friction which is very useful from the application point of 
view. From table 1 it is observed that the coefficient of skin friction increases due to increase in G„ P„ G c , co, 
K c , M and S whereas it decreases due to the presence of S c , K c and cot . 

From table 2 it is observed that the rate of heat transfer (N u ) at the surface increases in case of higher prandatl 
number (P r ), frequency of oscillation (co) and phase angle (cot) but it decreases due to the presence of heat 
source parameter (S). 

From table 3 it can be seen that the rate of mass transfer (S h ) at the surface increases with an increase in K c >0 
and S c but frequency (co) and phase angle (cot) has no significant effect. In case of generative reaction (Kc<0) it 
increases with an increase in the value of K c - 



Table - 2 Nusselt number (N u ) 
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Table - 3 Sherwood number (S h ) 
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IV. Conclusion 

Effect of elasticity on velocity depends on heating or cooling of the surface , being opposite to each other, 
whereas effect of other parameters are almost independent of heating or cooling of the plate. 

A flow characteristic in the present study is more dependent on mass buoyancy effect (G c ) rather than 
thermal buoyancy effect (G r ). 

Thinning of thermal boundary layer occurs for higher prandtl number fluid. 

Higher rate of mass transfer is experienced in case of heavier species and in the presence of destructive 
reaction. 
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